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Signatures of lensing of the cosmic microwave background radiation by gravitational potentials 
along the line of sight carry with them information on the matter distribution, neutrino masses, and 
dark energy properties. We examine the constraints that Planck, PolarBear, and CMBpol future 
data, including from the B-mode polarization or the lensing potential, will be able to place on these 
quantities. We simultaneously fit for neutrino mass and dark energy equation of state including time 
variation and early dark energy density, and compare the use of polarization power spectra with 
an optimal quadratic estimator of the lensing. Results are given as a function of systematics level 
from residual foreground contamination. A realistic CMBpol experiment can effectively constrain 
the sum of neutrino masses to within 0.05 eV and the fraction of early dark energy to 0.002. We also 
present a surprisingly simple prescription for calculating dark energy equation of state constraints 
in combination with supernova distances from JDEM. 



I. INTRODUCTION 

Precision studies of the cosmic microwave background 
(CMB) have helped us formulate a standard model of 
cosmology and measure several global parameters that 
describe our universe and its contents H, S 3- Six 
key parameters to describe the cosmology have been de- 
termined with 1-10% precision and CMB data plays a 
significant role in constraining other parameters, such as 
spatial curvature, the dark energy density, and the Hub- 
ble constant, in combination with other types of data. 

However, we know in some cases and allow the possibil- 
ity in other cases, that there are further fundamental pa- 
rameters beyond the six. One example is the mass of neu- 
trinos, where terrestrial experiments indicate a nonzero, 
though unknown jvalue: rrii, ^ 0.05 eV for at least one 
neutrino species ^ . Another set of parameters of great 
interest describes the properties of the dark energy caus- 
ing acceleration of the cosmic expansion. The dark en- 
ergy equation of state (EOS) may differ from the constant 
value 1/7 = —1 of the cosmological constant, and may vary 
with time. Indeed, this dynamics would be a key clue to 
the nature of the physics behind acceleration. The per- 
sistence of dark energy density to early times is another 
mystery that is crucial to explore. Current CMB data on 
the temperature and E-mode polarization spectra (and 
their cross-spectra) are of little use in themselves in ad- 
dressing these issues, and this holds to a large extent even 
in combination with other cosmological information such 
as supernova distances and large scale structure data. 

Fortunately, other types of CMB information exist, 
though they have not yet been measured. This includes 
the CMB deflection field - the action of gravitational po- 
tentials along the line of sight on the CMB - and the 
B-mode polarization spectra (and cross-spectra) result- 
ing from this. The effects carry contributions from all 
redshifts between the source (the last scattering surface 
at redshift z « 1090) and the observer, though like all 
lensing deflection the kernel from the geometric distance 
factors peaks approximately midway, z w 3 — 4. The 
growth of the gravitational potentials over this history 



carries within it information on the matter power spec- 
trum. Thus the effects of neutrino masses and dark en- 
ergy properties are encoded in the CMB. 

While subsets of these effects have been investigated 
before (see, e.g., [1, 0, i i [H), the effects have not 
generally been considered simultaneously (especially for 
dynamical dark energy) , with the critical covariances be- 
tween them. This article is also the flrst investigation of 
the important question of early dark density using CMB 
lensing. We also examine for a range of cases the added 
leverage of lensing information extraction through use of 
the optimal quadratic estimator which utilizes the unique 
non-Gaussian structure in the map caused by lensing. 

In S|TT]we lay out the methodology for obtaining preci- 
sion theoretical predictions for power spectra, and their 
slight variations with cosmology, and summarize the ob- 
servational capabilities of three benchmark CMB sur- 
veys. We explore adding neutrino mass to the standard, 
cosmological constant universe in mill and include as well 
the dark energy EOS and its time variation in fjTV] Dis- 
cussion includes complementarity with other cosmologi- 
cal probes and issues of foreground noise. In fjV] we in- 
vestigate early dark energy density, and present a simple 
prescription for cosmological constraints in WII We sum- 
marize the key prospects for intermediate range CMB 
experiments in Will 



II. POWER SPECTRA MODELING: THEORY 
AND EXPERIMENTS 

Primordial perturbations in the photon number den- 
sity arise from Gaussian, random, adiabatic fluctuations 
seeded in the inflationary era. These induce a photon 
temperature power spectrum, and interaction with infla- 
tionary gravitational waves and scattering from electrons 
creates _B-mode and i?-mode polarization power spectra 
(as well as a TE cross-spectrum), respectively. Grav- 
itational lensing shuffles the photon pattern on the sky 
^^'^ contributes to each of these spectra, as well as 
transforming some of the i?-modes into _B-modes, intro- 
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ducing a coupling between the two. Beyond these power 
spectra, lensing imprints non-Gaussianity into the CMB, 
and the CMB trispectrum encodes information about the 
deflection field power spectrum, or mapping of the pho- 
ton positions, itself [isj . 



A. Theory 

Accurate codes exist for computing each of these power 
spectra, at least for the standard cosmology. We uti- 
lize CMBeasy, which already implements several useful 
extensions to further cosmological parameters, includ- 
ing neutrino masses and several classes of dark energy 
[ij, [iBl ■ We have crosschecked results (for constant dark 
energy equation of state) with another code, CAMB, to 
ensure accuracy. Numerical stability is crucial, because 
several groups of cosmological parameters are highly de- 
generate and the differences between the power spectra 
for different cosmologies can be small, so numerical noise 
can distort the results. We carry out parameter estima- 
tion through Fisher matrix analysis. For the precision 
future data we consider, this should provide accurate con- 
straints. We check for convergence of the final results for 
various step sizes of the cosmological model differencing. 

The set of parameters considered includes the standard 
ones of primordial perturbation amplitude and power 
law index n, optical depth r, physical baryon density 
^Ibh^, cold dark matter density ^Ich^ and dark energy 
density Qde- The Hubble constant is a derived parameter 
/i^ = (rifc/i^ -t- ^lch'^)/{^ — fide) under the assumption of 
spatial flatness. The physical matter density is LUm = 
flbh^ + Vlch? . Since neutrinos are known to have mass 
and this influences the lensing and other power spectra, 
we always include as a parameter the physical neutrino 
energy density Vli,h^ or equivalently the sum of neutrino 
masses '^rriu — %A{^,yh?) eV. 

Since no guarantee exists that dark energy is a cosmo- 
logical constant, and generically other models have time 
variation of their equation of state, we consider two pa- 
rameters. Wo and Wa, to describe the dark energy equa- 
tion of state, w{a) = wq + u'a(l — a). Consideration of 
the physics behind dark energy led to this form and 
it has been shown to be accurate to 0.1% in describing 
observables [13] ■ The ACDM model corresponds to fixing 
Wo = —1, Wo = 0. Given that the CMB has strong sensi- 
tivity to the early universe, we also consider another class 
of dark energy models, early dark energy, where the dark 
energy density is non-negligible around and before the 
recombination epoch. These also have two parameters, 
the equation of state today wq and the constant high red- 
shift early dark energy density Vie [III- For z ^2 these 
look identical to the vjQ-Wa model where Wa ~ Sfig [3, 
but have distinct and possibly significant effects at high 
redshift. 

Thus we simultaneously fit either seven or 
nine parameters. We use the following fidu- 
cial parameter values throughout the paper: 



{A,,n,T,nbh'^,neh'^,nde:Y.^A = {2-41 X 
10-^ 0.963, 0.084, 0.02255, 0.1176, 0.72, 0.28 eV}. 



B. Deflection Field 

The angular power spectrum of the CMB has been 
used to constrain cosmological parameters with unprece- 
dented accuracy (see e.g. [l|), but its ability to inform 
us about the low redshift universe is limited by the so- 
called geometrical degeneracy. This arises because only 
angles are measured and, given some spectrum of primor- 
dial fluctuations, the power on each scale is nearly fixed 
for constant ^Jujm ^iss (where diss is the angular diame- 
ter distance to the CMB last scattering surface), which 
is degenerate under certain combinations of late universe 
parameters. (An exception to this arises on large angu- 
lar scales because the integrated Sachs- Wolfe (ISW) ef- 
fect [131 leaves another signature of dark energy on large 
scales, however owing to cosmic variance this effect is 
of limited use.) Also, the primordial CMB probes the 
baryon distribution at last scattering, which is smoothed 
on scales smaller than ~ 10' because of Silk damping 
f^H in the last scattering surface, while massive neutri- 
nos mostly impact matter agglomerations on projected 
smaller scales. 

The geometrical degeneracy can be broken by adding 
for example Type la supernova (SN) distance informa- 
tion or constraints on the expansion rate to the CMB 
power spectrum constraint (see, e.g., [l|). The effect of 
neutrinos on small scale structure can be pro bed through 
galaxy clustering or the Lyman-a forest j22l . |2^ . Alter- 
natively, or in addition, deflection of CMB photons on 
their way to us changes the statistics of the primordial 
pattern in a characteristic way that can be used to infer 
the lensing effect. What was originally a nearly Gaussian 
random field becomes non-Gaussian with the coherent 
correlation of patterns around large scale matter fluctu- 
ations. This type of non-Gaussianity, on a typical scale 
of approximately 2 degrees, is different from that used to 
study inflationary models [3, HI] , in that its three point 
function vanishes on most scales (except for those large 
scales on which the unlensed CMB is correlated with the 
lenses through the ISW effect). 

Lensing is described by the displacement vector of 
CMB photons on the sky, a(0), which is given as a{9) = 
£'cMB--Picn3 ^t0\ terms of the deflection angle 

d'^'^i^') i^r^p . S(x)^ j dDp(^,D), 

(1) 

where D is the angular diameter distance. The vector x 
describes the position in the lens plane, and the surface 
mass density (lensing can be imagined to good approx- 
imation as progressing through multiple, infinitely thin 
planes) is ^ projection of the three-dimensional 

density field p(x, D). In the so-called Limber approxima- 
tion, the lensing power spectrum C£" becomes a simple 
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integral over the matter power spectrum at all redshifts 
weighted by angular diameter distance ratios. In this pa- 
per we refer to modes in the lensing power spectrum as 
L and 1, to distinguish them from the CMB multipole I. 

Lensing also affects the angular power spectrum of the 
CMB The characteristic acoustic oscillation fea- 

tures are smeared out, as characteristically sized hot or 
cold spots are magnified or de-magnified by intervening 
lenses. The amount of over-smearing is scale dependent, 
encapsulating information about the shape of the mat- 
ter power spectrum, which in turn is affected by dark 
energy properties and neutrino masses. Because of the 
distance factors (the geometric kernel) and the growth 
factors, the matter power spectrum is best probed over 
the range z « 1 — 4. 

The effect of lensing on the CMB power spectrum is 
calculated within CMBeasy. In the presence of lensing, 
the power spectrum variance is not of the trivial Gaussian 
random field form. The non-Gaussian covariance is neg- 
ligible in temperature and E-mode polarization because 
the relative effect of lensing on these is small, however it is 
large for B-mode polarization, a factor of a few 0. The 
effects of marginalization when constraining individual 
parameters generally overwhelm the effect of the excess 
covariance however [2^ . We confirmed that the effect 
of non-Gaussian covariance on the parameter constraints 
in the next sections is negligible by checking that the un- 
certainties change by less than 10% (typically less than 
1%) if the sample variance in the B-mode is increased by 
a factor of five. 

The power spectrum over-smearing method provides 
a statistical estimate of lensing that is prone to sample 
variance because the actual distribution of the lenses on 
the sky remains unknown. To reconstruct the lensing po- 
tential ip (the line of sight projection of the gravitational 
potential of which the deflection vector a is the gradient) 
one needs to use the non-Gaussian information imprinted 
into the CMB. Lensing conserves surface brightness, so 
the probability distribution function of temperatures re- 
mains unchanged. Therefore the lowest order non-zero 
estimator of the lensing potential is quadratic. This has 
been investigated by |27l. 2811 and the minimum variance 
estimator was given by [29|. A quadratic estimator is 
generally of the form 



^(L)=iV(i) / 00(1)0'(L -1)5(1, L- 



1), (2) 



where 9 and 9' stand for temperature and/or polariza- 
tion modes on the sky, i.e. 9,9' = T,E,B. The optimal 
weight g and normalization N can be found using the fact 
that the deflected position can be written as a first or- 
der expansion of the displacement around the undeflected 
position, 6lL(x) = 6'UL(x-ha) = 6'UL(x) V^'(/'(x)V,6l(x). 
For the TT estimator, requiring an unbiased estimate 
and minimizing the variance leads to weighting of modes 



where Ci (Ci) is the unlensed (lensed) temperature power 
spectrum, following the convention of^ [3l|. Similar ex- 
pressions follow for polarization. The superscript "tot" 
originates from the fact that the lensed CMB and noise 
enter in the variance calculation. 

With the definition in Eq. ^ the noise of the lensing 
reconstruction equals the normalization which becomes 
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1 [(L - 1) • LCiL.ii + 1 • LG] 9(1, L - 1) 



(4) 

Physically the noise is a combination of instrumental and 
intrinsic shape noise (see below). 

Note that this is only the best quadratic estimator. 
Maximum likelihood methods can in principle be applied 
(S^ . [33} but they have been shown to only give small 
improvements for temperature and polarization experi- 
ments with the sensitivity levels assumed in this work, 
so we do not consider them here. We also note that the 
approximation above leads to a bias in the quadratic esti- 
mator, however for experiments considered in our paper, 
with angular resolutions larger than 3' as well as noise 
levels down to a micro-Kelvin, these are only a few per- 
cent and well understood (see [3^ 13^]). 

As is the case with the lensing of background galaxies, 
CMB lensing obtains most information from the small- 
est scale resolved by any given experiment as these al- 
low averaging over many background features. Because 
shapes in the CMB temperature can be intrinsically ellip- 
tical, averaging over many patterns becomes necessary to 
constrain relatively large lens features. Since unlensed B- 
type polarization patterns should be absent on scales less 
than a degree or so in concordance cosmology, quadratic 
estimators involving B, in particular the EB pair due 
to its higher signal-to-noise, are intrinsically more useful 
than temperature (as long as B can be imaged) and can 
be used to constrain lenses out to smaller scales. There- 
fore experiments beyond Planck, with the capability of 
imaging B-patterns, allow for reconstruction of lenses out 
to smaller angular scales [s^, HI] . 

In the following sections we will compare the lensed 
power spectra method (i.e. the over-smearing of acous- 
tic peaks) of inferring late universe parameter values to 
the optimal quadratic estimator (OQE) method. In the 
latter case we will use constraints on the unlensed power 
spectra in conjunction with a forecasted constraint on 
the lensing potential power spectrum^ Cf^ using Eq. (|4]) 
so we do not count the lensing information twice. 



5(1,L-1) = 



(L - 1) ■ LqL_i| -H ■ LC; 



(3) 



^ Note that other papers, for example l29l I30II . use the opposite 
notation to distinguish between lensed and unlensed spectra. 

^ Using the lensing potential power spectrum is equivalent to us- 
ing the deflection power spectrum. They are simply related by 
= L^cf^^ (in the flat sky approximation applied here). 
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C. Experiments 

In this paper we consider three different experiments, 
two of which are scheduled to begin observations in the 
near future, to forecast constraints on neutrinos and dark 
energy. The Planck satelhte will be launched in the 
second quarter of 2009 and will observe the full sky from 
the semi-stable Lagrange point L2. We take into account 
a foreground cut for galactic emission and assume a sky 
coverage of 0.75 to be useful for cosmological analysis. 
We have adopted the experimental specification values 
in @. 

Combining both large sky coverage and high sensitiv- 
ity, we consider the futuristic CMBpol concept of a 
satellite mission specialized on polarization with ultra- 
high sensitivity. We have used values from Our 
assumed specifications are summarized in Table ID We 
postpone further discussion of PolarBear, an interme- 
diate term and sensitivity experiment until Section [VIII 



Experiment 


V 


/sky 


SpWHM 


At 


Ap 


Planck 


100 GHz 


0.75 


9.2' 


51 






142 GHz 


0.75 


7.1' 


43 


78 




217 GHz 


0.75 


5.0' 


65 


135 


PolarBear 


150 GHz 


0.025 


4.0' 


3.5 


5 




220 GHz 


0.025 


2.7' 


8.5 


12 


CMBpol 


all freq. comb. 


0.75 


3' 


1 





TABLE L Experimental specifications assumed in the fore- 
casts in this paper, for the various frequency bands of 
Planck, PolarBear, and CMBpol. The temperature and 
polarization sensitivities At, Ap are given in units of /xK- 
arcmin. 

From these experimental characteristics the full esti- 
mator covariance matrices for each multipole / can be 
constructed (e.g. jlll). The (Gaussian) covariances be- 
tween the power spectrum and cross correlation estima- 
tors are given by 

(C?'"' + nP'') + (C,^'" + N™) {CY^ + nY^)] (5) 

where the noise power spectrum'^ 

for XX = TT, EE, BB, and xf''' is given by Eq. (g]) 
and X^^ = when X ^Y. Here At and = As = 



^ When there are multiple frequency bands, the total noise power 
spectrum is given by N^^^^ = where the sum is over 

the individual bands and we have suppressed the superscripts. 



Ap are the temperature and polarization sensitivities, 
^FWHM is the angular resolution, and Tq is the tempera- 
ture of the CMB today. 

Figure [T] shows that the noise of the experiments con- 
sidered here is so low compared to the signal that they 
gather much of their information from scales beyond 
I — 2000 in the temperature power spectrum (Polar- 
Bear curves, not shown, would lie between Planck and 
CMBpol curves). This is especially true for lensing, be- 
cause the characteristic displacement of a CMB photon 
on its way from the last scattering surface to us is of 
order 2-3 arcminutes, and the smallest scale resolved by 
a given experiment contains most of the lensing infor- 
mation. However on scales I > 2000 in the temperature 
power spectrum secondary anisotropics that are larger in 
magnitude than lensing, such as the Sunyaev-Zel'dovich 
(SZ) effects [S^l and radio as well as infrared point sources 
will make extraction of lensing information challenging. 
This is true as well for the optimal quadratic estimator, 
which might get confused by the extra non-Gaussianity 
carried by these foregrounds. In addition to the limi- 
tation due to the instrumental noise level and angular 
resolution, we therefore also quote our results with high- 
l cuts at different scales, to show how these foregrounds 
affect parameter constraints. We note that while point 
sources and the SZ are expected to be significantly dim- 
mer in polarization than in temperature [381 . |3S] , there 
the cutoff at high I does not lead to as much loss in infor- 
mation as the polarization signal-to-noise ratio is small 
on angular scales beyond I — 2000. 

Finally, the Fisher matrix is given by the expectation 
value of the second derivative of the logarithm of the 
likelihood function C{Ci\9i). Assuming Gaussianity of 
the likelihood it is of the form 

I a.J3 

where a and P run over the five observables: temper- 
ature, E-mode polarization, T-E cross correlation, B- 
mode polarization, and lensing potential power spectrum 
(where the OQE is used), and i,j run over the cosmolog- 
ical parameters. The covariance matrix between param- 
eters is given by the inverse of the Fisher matrix. 



III. NEUTRINO MASS CONSTRAINTS IN 
ACDM 

We begin looking at cosmological constraints in the 
simplest model consistent with both cosmological and lo- 
cal observations: a cosmological constant universe with 
non-zero mass neutrinos. Three types of data cuts are 
employed - by classes of observations, experiments, and 
systematics. 

The classes of observations are 1) unlensed TT, TE, 
EE power spectra, 2) adding the effect of lensing to 1), 
3) adding the BB power spectrum to 2), and 4) using 
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FIG. 1; Temperature and polarization power spectra TqI{1 + l)Ci/2iT[fiK^] vs. multipole / for the ACDM fiducial cosmology. 
TT, EE, BB, and TE spectra (solid curves, dotted where negative) run clockwise from upper left. Dashed curves show the 
power spectrum errors, TqI{1 + l)ACi/2-K[jj.K'^], for the Planck (long dash) and CMBpol (short dash) experiments. 



1) plus information on the lensing potential through the 
optimal quadratic estimator discussed in mi Bl This al- 
lows understanding of the effects of lensing on just the 
temperature and E-mode spectra, the information in just 
the BB power spectrum caused by lensing, and methods 
for using the complete effects of lensing. 

On the experimental side, we consider Planck, slated 
for launch in mid-2009, and the far future CMBpol mis- 
sion. Discussion of the impact of intermediate scale 



ground-based missions is postponed until Section IVIII 
Additionally we examine the influence of the level of sys- 
tematics in terms of ^max, such as induced through fore- 
grounds external to the experiments. 

Figures [2] and [3] illustrate the constraints on neutrino 
mass and dark energy density (cosmological constant) for 
the different data set types and systematics levels. All 
figures show 68% confidence level contours; the fiducial 
model is ACDM, with ^m^, = 0.28 eV. Use of lensing 
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information clearly adds substantial leverage, and mea- 
surement of B-modes or the lensing potential play an 
important role. The two methods of including the full 
lensing information - adding B-modes or adding the lens- 
ing potential - are nearly equivalent (see Will for further 
discussion of this). 

Considering the constraints for different systematics 
levels, we see that much of the lensing leverage is achieved 
by ^max ~ 2000. On smaller scales point sources and 
the SZ effects are expected to dominate over lensing and 
our limited ability to clean foregrounds through multi- 
frequency observations will likely not allow lensing re- 
construction on much smaller scales. 

Finally, the dramatic improvement of CMBpol over 
Planck is clear in Fig. 01 Here we adopt as a standard 
systematics limit Zmax = 2000 and show the confidence 
contours for each data set type for both experiments. 
While lensing information does improve the Planck con- 
straints, it runs into a wall due to the relatively high 
instrumental noise. Furthermore, Planck essentially can- 
not see B-mode lensing at all (see Fig. [T]). This is one 
of the motivations for intermediate experiments such as 
PolarBear. 



IV. ADDING DARK ENERGY DYNAMICS 

The cosmic microwave background plays a crucial role 
in breaking the degeneracies of other probes in order to 
constrain the properties of dark energy. However, un- 
lensed CMB data itself has very little leverage on learning 
about dark energy, since the power spectra reflect mostly 
conditions in the high redshift universe or at best a sin- 
gle weighted average of dark energy influence through the 
distance to last scattering. With the addition of lensed 
CMB data we can ask if this improves the leverage on 
dark energy; we emphasize that it is crucial to consider 
at least minimally realistic models that include dynam- 
ics in the EOS: taking the value of w constant from the 
present to z « 1100 is highly non-generic. 

It is also important to retain the inclusion of neutrino 
mass while making this investigation; both neutrino mass 
and dark energy influence the CMB in many of the same 
ways, e.g. suppressing structure and causing gravitational 
potentials to decay. Ignoring neutrino mass could lead to 
overoptimistic constraints on dark energy. In this section 
therefore we add wq and Wa as fit parameters to the set 
considered in the previous section. We explore the con- 
straints under the same variety of data cuts as in that 
section. 

The geometric degeneracy due to the acoustic peaks 
feeling dark energy mostly through the integrated dis- 
tance to last scattering remains strong, and no reason- 
able constraints can be placed on the dark energy EOS 
even with full use of the lensing information. We there- 
fore turn to the issue of complementarity: does the CMB 
data substantially help other probes of dark energy? In 
particular we examine complementarity with luminosity 



distances measured by Type la supernovae, since the two 
probes are well known to strengthen each other [13, EH . 
We consider luminosity distances measured to ~ 1% from 
z — — 1.7, including systematics, as could be provided 
by a supernova sample realized by a SNAP-type Joint 
Dark Energy Mission (JDEM)[4l|. 

Figure O shows the constraints in the wo-Wa plane, 
marginalizing over the other seven parameters, for each 
data set type. The first thing to notice is the clear im- 
provement in measuring the time variation Wa over the 
supernova sample alone due to even unlensed CMB data. 
Adding lensed CMB data continues to tighten the con- 
straints, in both Wo and Wa, except for the worst system- 
atics level Zinax = 1000. Full lensing information contin- 
ues the improvement modestly on the limits, and some- 
what narrows the contours. 

Figure [6] exhibits the analogous situation for different 
systematic limits ^max- In contrast to the ACDM case, 
here the constraints continue to improve for higher Zmax- 
There are also slight differences between the two methods 
of fully incorporating lensing: use of B-modes or OQE of 
the lensing potential. This emphasizes that conclusions 
on systematics or analysis methods should not be based 
solely on examination of the vanilla ACDM cosmology. 
Finally, when systematics are low, Zmax — 4000, sufficient 
information is present in the lensed E-modes that further 
lensing information is unimportant. 

The improvements in dark energy estimation that 
CMB lensing brings is illustrated in Fig. [7] as a function 
of experiment. For Planck, again no lensing informa- 
tion beyond E-modes is useful, though the contour area 
decreases by a factor 1.9 from the unlensed case to the 
OQE case. By contrast, CMBpol could reduce the likeli- 
hood contour area by a factor 4.2 relative to the unlensed 
Planck case, with the full lensing information helping by 
a factor 2.7 relative to unlensed CMBpol. 

To test the effect of including both neutrino mass and 
dark energy dynamics. Fig. [5] shows the likelihood con- 
tours for the Zinax = 2000 CMBpol case, marginalizing 
over vs. fixing ^ rrii, . The fully marginalized uncertain- 
ties are a{Y,m,,) = 0.041, a{wo) = 0.066, a{wa) = 0.25. 
While the la limits on the parameters do not change that 
strongly, the total area of the contour is significantly af- 
fected. For the unlensed (fully lensed) case the area in- 
creases by a factor 2.9 (1.5) when properly marginalizing 
over neutrino mass. (This effect would be more severe 
when considering CMB data alone.) Note that for the 
CMBpol case the correlation coefficient between wq and 
is 0.23 and between Wa and J^''^'^ —0.41; while 
not highly correlated, these are sufficient to give the ap- 
preciable effect. 



V. EXPLORING EARLY DARK ENERGY 

In ACDM, the fractional contribution of dark energy 
density is of order 10~^ at last scattering. However, many 
models exist where this can be at the percent level j43| . 
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FIG. 2: Cosmological constraints on the neutrino mass and dark energy density in the ACDM fiducial cosmology from CMBpol. 
Within each panel the contours correspond to systematic cuts at Zmax = 1000, 2000, 3000, 4000 from outer to inner. The panels 
use different data cuts: no lensing (upper left), including lensing from T- and E-modes (upper right), including lensing from 
T-, E- and B-modes (lower left), and including lensing through the optimal quadratic estimator of the lensing potential (lower 
right). 



vifith important impacts on the sound horizon scale and 
baryon acoustic osciUations, structure formation, and 
secondary anisotropics [H, [3, [H, HI, [Ulii, . Such 
early dark energy models follow from physics where the 
dark energy traces the energy density of the dominant 
component of the universe, as in high energy physics and 
string theory models with dilatation symmetries j48| . 



Although the sound horizon is altered in the presence 
of early dark energy by ^ (1 — Q^Y^^, this shift can 
be hidden in the temperature power spectrum by com- 
pensating changes in the other parameters [l^. This is 
problematic for baryon acoustic oscillation experiments, 
which use the sound horizon as a standard ruler to probe 
cosmology through distances. More generally, definitive 
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FIG. 3: As Fig.[2]but here within each panel the contours correspond to data set types, and the panels use different systematics 
levels: Z^ax = 1000 (upper left), 2000 (upper right), 3000 (lower left), 4000 (lower right). 



recognition of early dark energy is quite important to 
have confidence in the accurate estimation of the other 
parameters, ensuring that they are not biased due to in- 
correctly assuming no early dark energy. Furthermore, 
detection of early dark energy would immediately give 
crucial clues to understanding the nature of dark energy. 

Since CMB lensing depends on the growth of struc- 
ture, it is a good candidate for constraining dark energy. 
More generally, hints already exist in [l9| that polariza- 
tion information can help break degeneracies involving 



early dark energy. Here we carry out a more comprehen- 
sive likelihood analysis for unlensed polarization power 
spectra and examine for the first time CMB lensing con- 
straints on early dark energy. To do this, we employ the 
parametrization for the fractional dark energy density as 
a function of scale factor proposed by [ll]. 



^de — J^e (1 



-3uio'> 



rip (l — a 



(8) 



where fide is the current dark energy density, is the 
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trino mass and dark energy density in the ACDM fiducial 
cosmology from Planck (dashed contours) vs. CMBpol (solid), 
taking Zmax = 2000. 



constant dark energy density at early times, and wq is the 
present dark energy equation of state. Hence, the two 
added parameters Qe and wq describe the dark energy 
properties. 

Figure [5] shows the constraints in the wo-^^e plane, 
marginalizing over the other seven parameters, for differ- 
ent data set types. The fiducial model has wq — —0.95, 
He = 0.03. As in the WQ-Wa case, the CMB degeneracies 
are too strong to allow constraints by the CMB alone, so 
we have again folded in supernova distance data (which 
does not directly constrain fig). We see that unlensed 
power spectra including polarization information can in- 
deed tightly constrain early dark energy. Adding lensed 
CMB information in fact mostly constrains further wq, 
having minimal effect on ^lf,. Recall from i jll Al that out 
to 2; w 2, the early dark energy model looks very much 
like a standard WQ-Wa model that would not give appre- 
ciable early dark energy density. Thus, early dark energy 
is too early for even the broad redshift kernel of CMB 
lensing to have significant sensitivity to it. 

Figure [TO] exhibits the analogous situation for differ- 
ent systematic limits /max- Again in contrast to the 
ACDM case, here the constraints continue to improve 
for higher Z^ax, although less rapidly for Z^ax ^ 3000. 
The fully marginalized uncertainties for the Zmax = 2000, 
full lensing case are (t(X]"t-i/) = 0.047, cr{wo) — 0.018, 
a{ile) — 0.0019. This is an impressive constraint on the 
early dark energy density, able to give definite guidance 
to the nature of dark energy, ruling out classes of models. 



Because CMBpol would have much better polarization 
measurements than Planck, it will constrain f2e better 
by a factor 2.2, as shown in Fig. [TT] The area of the 
dark energy properties' confidence contour improves by 
a factor 3.9. 

Finally, we summarize our results for the dark energy 
and neutrino mass uncertainties in Table ITTl for the three 
cosmological models considered, assuming Zmax — 2000. 
However, one should see the figures for the full contours. 
Due to degeneracies in the presence of dynamical dark 
energy, we add supernova data in these cases to con- 
strain the dark energy equation of state, although the 
uncertainties on fig and ^ are not strongly affected. 



Model Experiment 


a{wo) 


0-{Wa) 


(j(ne) 


(T(Em„) [eV] 


ACDM Planck 








0.11 


ACDM CMBpol 








0.036 


WQ-Wa Planck-|-SN 


0.073 


0.32 




0.13 


wo-Wa CMBpol-fSN 


0.066 


0.25 




0.041 


wo-fie Planck-I-SN 


0.032 




0.0041 


0.15 


wo-fie CMBpol-FSN 


0.018 




0.0019 


0.047 



TABLE 11: Uncertainties in parameters beyond standard 
ACDM for Planck and CMBpol. In all cases, we use un- 
lensed temperature and polarization spectra and the optimal 
quadratic estimator of the lensing spectrum to extract cosmo- 
logical information from the CMB data. For cases involving 
dynamical dark energy we fold in supernova distance infor- 
mation from a SNAP-like JDEM experiment, although this 
mostly affects only the uncertainties on wo, Wa- 



VI. SHORTCUT FOR JOINT DARK ENERGY 
CONSTRAINTS 

As seen in Section HVl when CMB and supernova data 
are combined, we can obtain strong constraints on the 
nature of dark energy. While the supernova data depen- 
dence on cosmological parameters is straightforward, cal- 
culating a CMB Fisher matrix can be quite time consum- 
ing. The procedure requires computing multiple CMB 
spectra using a Boltzmann code (CMBeasy in our case) 
for different values within a set of cosmological parame- 
ters in order to obtain the derivatives of the observables 
with respect to the cosmological parameters. 

To investigate a range of cosmological models it would 
therefore be quite useful to have a shortcut to calcu- 
lating the constraints on the dark energy parameters 
WQ,Wa,HDE from CMB data. One such shortcut is his- 
torically well known, the shift parameter A9l to encap- 
sulate the information in the temperature power spec- 
trum acoustic peaks. However, as polarization data gets 
added, other parameters have been suggested as addi- 
tions, e.g. the acoustic peak scale I a [H, although (soj 
showed that the shift parameter is still quite accurate. 
Here we investigate the cosmological constraints from 
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FIG. 5: Cosmological constraints on the dark energy equation of state parameters wo and Wa from CMBpol in combination 
with SNAP-quality supernova distances. Within each panel the contours correspond to systematic cuts at /max ~ 1000, 2000, 
3000, 4000 from outer to inner. The panels use different data cuts: no lensing (upper left), including lensing from T- and 
E-modes (upper right), including lensing from T-, E- and B-modes (lower left), and including lensing through the optimal 
quadratic estimator of the lensing potential (lower right). The dotted curve gives the constraints from supernovae alone. 



combining CMB temperature, polarization, and possibly 

deflection, spectra and supernova data, and we show that 
a simple use of the shift parameter has excellent acciiracy. 

Specifically, for constraints on the dark energy param- 
eters a strong prior on the shift parameter, or reduced 

distance to last scattering, d = ^uj,n c^iss j is nearly equiv- 
alent to the full CMB data, even including polarization 



and lensing data. That is, the CMB Fisher matrix for 

nDE,W(),Wa after marginalizing over the other param- 
eters is almost identical to the Fisher matrix calculated 
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FIG. 6: As Fig.[5]but here within each panel the contours correspond to data set types, and the panels use different systematics 
levels: imax = 1000 (upper left), 2000 (upper right), 3000 (lower left), 4000 (lower right). Since using lensed TT/EE/TE spectra 
is not a matter of simply adding to the Fisher matrix from unlensed spectra, it is possible for a lensed contour to lie slightly 
outside of the unlensed contour, as in the Zmax ~ 3000 case. 



from a single constraint on d ^ . The prior on the quantity 
d required to match the CMB data depends on the CMB 
experiment and on whether or not we fix the neutrino 



mass. We emphasize that the level of the prior does not 
correspond to the actual determination of d from the ex- 
periment, because the prior also encodes the other spec- 
tra information. For the CMB experiments we consider, 
the equivalent prior on d is 0.2% — 1.2%. 



Note this holds for the CMB Fisher matrix itself, without any 
supernova information. 



Note that because early dark energy does not merely 
affect the projection of the last scattering surface onto 
our sky, but also affects the shape of the anisotropy spec- 
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FIG. 8: As Fig. [7]for CMBpol only, but here showing the ef- 
fect of fixing (dashed contours) rather than marginaliz- 
ing over it (solid) as is standard for all parameters not shown. 



trum at last scattering directly, we do not expect the d 



prior to be a complete description there and indeed we 
found the prior is not effective in this case. 

We compare the shift parameter prescription to the use 
of the actual CMB Fisher matrix in Fig. [12] by consider- 
ing la joint contours in the wq ~ Wa plane for CMB -t- 
SN. For Planck (top panel), if we marginalize over ^ 
and if we do not include the information from CMB lens- 
ing, the constraints from the CMB-fSN are almost ex- 
actly the same as those with a 1.2% prior on d. The 
constraints are improved quite a bit if the lensing infor- 
mation is added. In this case, the constraints are about 
the same as the constraints one gets with a 0.6% prior 
on d. In the case of fixing ^ m,y instead of marginalizing 
over it, the shift parameter prior applies as well, at 0.2% 
matching the CMB-f-SN contours whether lensing infor- 
mation is used or not. Interestingly, in the case of fixed 
^ , adding lensing information does not appreciably 
improve the constraints on wq and Wa any more. 

Note that the Planck experiment can be approximated 
extremely well by the shift parameter prior in all these 
cases. The extent a{wo), a(wa), width a{wp), area 
1/VdetF, and orientation of the dark energy EOS con- 
tours match, as seen in Fig.[n]and quantified in TableHTTl 



Data 


a{wo) 




a{wp) 


VdetF 


SN-KPlanck 


0.073 


0.32 


0.031 


101 


SN-f 0.6% d 


0.076 


0.31 


0.032 


99 


SN-hCMBpol 


0.066 


0.25 


0.018 


223 


SN-hO.2% d 


0.076 


0.29 


0.017 


202 



TABLE III: Dark energy constraints from supernovae and 
CMB compared to constraints from supernovae and a prior 
on the shift parameter d. We assume imax = 2000 and use 
the optimal quadratic estimator to extract lensing informa- 
tion for both Planck and CMBpol. We marginalize over the 
sum of the neutrino masses and over the other parameters 
of the model. Note a{wp) is the width of the WQ-Wa contour 
at Wa = (i.e. the uncertainty in constant w) and VdetF is 
the inverse area of the contour (sometimes used as a figure of 
merit). 



For CMBpol (bottom panel of Fig. fT^ . the constraints 
on dark energy can be very well approximated by a 0.2% 
prior on the shift parameter. This is true independent of 
whether one fixes or marginalizes over it because for 
CMBpol with lensing, fixing only improves the con- 
straints on wq and Wa a little bit compared to marginal- 
izing over rrii, (see Fig. |8]) . Note that making the prior 
on d even smaller than 0.2% does not change the contour 
significantly. To illustrate this. Fig. [T2l shows the contour 
for a prior of 0.0001%, essentially fixing d. It is almost 
the same as the contour for 0.2%. 

The combination of supernova data with a prior on d 
always gives an ellipse with ends touching the contour 
from supernovae alone. This means that while both the 
area enclosed by the contour and the uncertainty in Wa 
may be improved, the uncertainty in wq is the same as 
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FIG. 9: Cosmological constraints on the early dark energy fraction ile and present equation of state parameter wo from 
CMBpol in combination with SNAP-quahty supernova distances. Within each panel the contours correspond to systematic 
cuts at ^max ~ 1000, 2000, 3000, 4000 from outer to inner. The panels use different data cuts: no lensing (upper left), including 
lensing from T- and E-modes (upper right), including lensing from T-, E- and B- modes (lower left), and including lensing 
through the optimal quadratic estimator of the lensing potential (lower right) . 



the uncertainty from supernova data only. Since Planck 
constraints are described almost perfectly by the shift pa- 
rameter, this is also true for Planck. However, once we 
include precision measurements of polarization by con- 
sidering CMBpol, the ends of the error ellipse can move 
away from the "SN only" contour and thus (slightly) im- 
prove the constraint on wq. This effect cannot be re- 



produced by the prior on the shift parameter. Hence, 
for CMBpol, the shift parameter prescription works less 
well than for Planck, although it is still quite adequate. 
Again, Table [HT] quantifies the accuracy of substituting 
the prior in place of the full CMB spectra. 
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FIG. 10: As Fig.[9]but here within each panel the contours correspond to data set types, and the panels use different systematics 
levels: Zmax = 1000 (upper left), 2000 (upper right), 3000 (lower left), 4000 (lower right). Since using lensed TT/EE/TE spectra 
is not a matter of simply adding to the Fisher matrix from unlensed spectra, it is possible for a lensed contour to lie slightly 
outside of the unlensed contour. 



VII. PROGRESS IN NEAR-TERM 
EXPERIMENTS: POLARBEAR 



In this section we explore the merit of near term 
ground-based polarization sensitive CMB missions to 
constrain dark energy and neutrino properties. A num- 
ber of such experiments are currently being built or have 
been funded including B1CEP/BICEP2 BRAIN 



[12, CfOVER t53i], EBEX [H, QUIET [5|, Spider [s^, 
SPTpol. Here we focus on one of them, PolarBear, 
as it represents a good combination of the high angular 
resolution and sensitivity some of these experiments will 
be capable of. 

PolarBear is a ground based telescope with sched- 
uled beginning of operations in 2009, and deployment 
to Northern Chile in 2010. It plans to observe 2.5% of 
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FIG. 11: Comparing the cosmological constraints on the early 
dark energy fraction Qe and present equation of state param- 
eter Wo from Planck (dashed contours) vs. CMBpol (solid), 
taking Zmax = 2000 and including SNAP-quality supernova 
distances. 



the sky. The low noise of its detectors will enable this 
experiment to go beyond Planck in imaging the B-type 
polarization pattern, which on small scales is a clear sig- 
nature of gravitational lensing as it cannot be produced 
by scalar fluctuations. However, the smaller sky cover- 
age does not allow the lensing potential power spectrum 
to be constrained with as high a signal-to-noise on most 
scales, making forecasted constraints generally somewhat 
less good. To describe this experiment's capabilities, we 
have adopted specifications from [H^], and the resulting 
likelihood contours are shown in Fig. [T31 

The cosmological constraints from PolarBear lensing 
reconstruction are less good than those from Planck, de- 
spite the significantly lower noise level. The reason is 
simply that the limited sky coverage does not allow most 
modes in the temperature, polarization, and lensing po- 
tential power spectra to be constrained with as high over- 
all signal-to-noise. However the constraints are still in- 
teresting relative to current limits. Moreover, we partic- 
ularly note that our parameter space has been limited 
to not include tensor fluctuations, which are a natural 
consequence of inflationary models. With its low noise 
level PolarBear will attempt to measure these gravita- 
tional waves from inflation and will help break degenera- 
cies between the tensor-to-scalar ratio and other param- 
eters that are present in the Planck data. Furthermore, 
we have not included running of the scalar spectral index; 
again, PolarBear's high resolution and low noise will pro- 
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FIG. 12: Joint constraints from CMB and supernovae can be 
well approximated by simply replacing the CMB data by an 
appropriately sized prior on the shift parameter d — rfiss . 
[Top panel] Combining Planck data with supernovae, the dark 
blue (light red) curves represent not using (using) lensing in- 
formation. These two cases are well approximated by replac- 
ing CMB data by d priors of 1.2% and 0.6% respectively. 
[Bottom panel] Combining CMBpol data with supernovae, 
the light red curve represents usin^ lensing information. This 
is fairly well approximated by a d prior of 0.2% (outer black 
curve). Tighter priors have little effect (see inner black curve). 
Both panels take Zmax = 2000, and the dotted line in both 
panels is the contour from just supernova data. 
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FIG. 13: Constraints from the intermediate experiment PolarBear are not as strong as Planck within the restricted inflationary 
scenario assumed here. Contours are constructed using unlensed TT/TE/EE data plus the optimal quadratic estimator for the 
lensing spectrum. Blue dotted contours repeat the results for Planck from Figs. l4l [71 1111 



vide an advantage in breaking degeneracies once running 
is included. 

We have found that with Planck the use of the 
quadratic estimator vs. lensed power spectra leads to a 
significant improvement of the constraints on parame- 
ters to which lensing is sensitive. To be specific, we find 
a 39% improvement on the neutrino mass scale and a 
26% improvement on J^a- The improvement in the case 
of PolarBear and CMBpol is however only marginal. 



To illuminate this trend, in Figure [HI we plot the power 
spectra of the lensing potential and lensing reconstruc- 
tion noises as well as the total errors. The dotted lines 
show the lensing reconstruction noises for each experi- 
ment. PolarBear has better capability to map the lensing 
potential in the observed patches on the sky than Planck 
(although it reconstructs far fewer of these patches and 
therefore the total error is larger than for Planck). The 
lower lensing noise feeds into the estimation with the op- 
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FIG. 14: The lensing potential power spectrum (solid line) 
is shown together with the lensing reconstruction noises for 
the three experiments considered in this work (dotted lines) 
and the total error on individual multipoles in the lensing 
potential, a combination of sample variance and noise (dashed 
lines) . 



VIII. CONCLUSIONS 

Continued advances are expected in measuring the 
cosmic microwave background radiation including lower 
noise and better systematics control, smaller beams and 
wider surveys, and extension to polarization, cross spec- 
tra, and CMB lensing information. These will greatly 
improve our knowledge of a variety of cosmological pa- 
rameters related to primordial perturbations. Here we 
have explored their impact on physics where the CMB 
has not had as much direct leverage - extensions to the 
standard model of cosmology such as the necessary neu- 
trino mass and the suspected dynamics of dark energy. 

We find the following general points to guide the design 
and analysis of CMB experiments, both ground based 
and the CMBpol satellite concept: 

• Systematics, such as point sources and other fore- 
ground contamination, will affect the lensing po- 
tential and other power spectra, and should be re- 
moved at the level of at least Zmax = 2000. Con- 
straints improve only slowly for higher ^max when 
using the full information in the CMB. 



• Analysis of gravitational lensing of the CMB can 
proceed either through consideration of induced B- 
mode polarization or through an optimal quadratic 
estimator directly of the deflection field; the opti- 
mum is not steep so the two approaches are nearly 
equivalent for these purposes with data beyond 
Planck. 

• For exploration of suites of cosmological models, we 
establish the accuracy of a shortcut in terms of an 
effective prior on the CMB shift parameter. This 
is remarkably efficient in summarizing the informa- 
tion from the CMB spectra. 

Determination of the sum of neutrino masses can be ac- 
complished by CMBpol with an uncertainty of 0.05 eV, 
marginalizing over all other parameters including dark 
energy properties. This corresponds to greater than a Str 
detection for the fiducial value adopted, and represents 
a factor 3 improvement over Planck expectations. Re- 
stricted to a ACDM cosmology, the constraints tighten 
by a factor ~ 1.3. 

Determination of the dynamical properties of dark en- 
ergy is less powerful. Complementary information, such 
as from distance measurements, is required with the 
leverage of the two data sets together allowing significant 
constraints. The present dark energy equation of state 
Wo could be estimated to 0.07 and the time variation 
Wa to 0.25, including marginalization over other cosmo- 
logical parameters including neutrino mass. This would 
improve further as other probes are added. While the 
marginalized constraints do not improve greatly in go- 
ing from Planck to CMBpol, the area of the uncertainty 
contour shrinks by a factor 2. 

The most significant impact from the CMB comes 
within early dark energy models. Here the improvement 
from Planck to CMBpol is a factor 2 in estimation of 
both wq (to 0.02 for CMBpol plus distances) and early 
dark energy density fie (to 0.002 for CMBpol plus dis- 
tances), while the uncertainty area shrinks by a factor 4. 
This provides the possibility of a ~ lOtr detection of early 
dark energy, which would immediately revolutionize our 
physics thinking. 

CMB lensing offers an intriguing new window on the 
universe, especially because of its sensitivity to the prop- 
erties of expansion and growth in the poorly probed 
epoch z « 1 — 4. Experiments nearly in the process 
of data collection will teach us not only about the pri- 
mordial conditions but also about the interesting period 
when dark energy first becomes significant, as well as es- 
tablishing a link to terrestrial experiments to measure 
the neutrino masses. 
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